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Within the mean field theory we extend the eflective quasi-one-dimensional nonpolynomial 
Schrodinger equation (NPSE) approach to the description of a spin-1 atomic condensate in a tight 
radial confinement geometry for both weak and strong atom-atom interactions. Detailed compar- 
isons with full time dependent 3D numerical simulations show excellent agreement as in the case of 
a single component scalar condensate, demonstrating our result as an efficient and effective tool for 
the understanding of spin-1 condensate dynamics observed in several recent experiments. 
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Although our ability to perform numerical simulations 
keeps increasing with computer technology, full 3D time 
dependent calculations still represent a significant chal- 
lenge. In many situations, one explores the inherent sys- 
tem symmetries, e.g., cylindrical and spherical symme- 
tries in space, to reduce the number of spatial dimension 
from 3D to 2D or even ID. The description of atomic 
condensate dynamics in terms of a mean field theory is 
such an example. With a tight radial confinement, a con- 
densate becomes cigar-shaped. Several effective ID ap- 
proaches have been developed P, 0, Q , with the simplest 
of them assuming a fixed transverse Gaussian profile. Re- 
cent studies, however, have indicated that the effective 
quasi-one-dimensional (ID) nonpolynomial Schrodinger 
equation (NPSE) is the most powerful and efficient tool, 
at least for a weakly interacting atomic condensate [ij . In 
this brief report, we generalize such an NPSE approach to 
the case of a spin-1 atomic condensate in a cigar-shaped 
trap. 

Multi-component atomic condensates or spinor con- 
densates have become an actively investigated topic in 
atomic quantum gases |lllllllllli|litil|ll|li|ll 
lig . iwL [13 • The examples that have been experimentally 
realized include that of two-component pseudo-spin-1/2 
0, three-coniponent spin-1, and five-component spin-2 
condensates H, IE IE Ol- Several recent experiments 
have observed interesting coherent spatial fragmentation 
of the spin-1 condensate when it is confined in a single 
running wave optical trap, i.e., in a cigar-shaped trap. To 
provide a proper theoretical description for these obser- 
vations, numerical approaches have been used to study 
the nonlinear spatial-temporal dynamics for a spin-1 con- 
densate. It is therefore desirable to have a more efficient 
theoretical approach instead of the 3D coupled Gross- 
Pitaevskii (GP) equations that is uniformly valid to both 
strongly and weakly interacting limits. 

A spin-1 Bose condensate is described by the Hamil- 



tonian in second quantized form (repeated indices are 
summed) as 



n 



df 



2M 



V + Fext + ) (1) 
I , 



where ^^(t^) (*]) is the field operator that annihilates 
(creates) an atom in the jth internal state at location r, 
j = +, 0, — denotes atomic hyperfine state \F — l,mp — 
-f 1, 0, — 1), respectively. M is the mass of each atom and 
Vcxtir) is an internal-state- independent trap potential. 
Terms with coefficients cq and C2 of Eq. ^ describe 
elastic collisions of two spin-1 atoms, expressed in terms 
of the scattering lengths ao (02) in the combined symmet- 
ric channel of total spin (2), Co = 4:TTn'^{ao + 2a2)/3M 
andc2 = 47r?i^(a2-ao)/3M. =x,y,z are spin-1 matrices 
[TtI . Since an external magnetic field is usually present 
in experiments, it is also included in our formulation. For 
simplicity, the magnetic field B is taken along the quan- 
tization axis (z). The Zeeman shift on each atomic state 
is then given by the Breit-Rabi formula 0| . 

Adopting the mean field theory when the condensate 
consists of a large number of atoms, we introduce the 
condensate order parameter or wave function $i = 
for the ith component. Neglecting quantum fluctuations 
we arrive at the mean field energy functional, 
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from which the coupled GP equations can be derived ac- 
cording to ihd^i/dt = 5£/d^*. They are given below in 
explicit form as 
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where n — rii is the total condensate density and Ui = 

The external trap is assumed harmonic V^oxt — M{ijj\r\ + 101x^)12 with cylindrical symmetry ujx — ujy — uj±, and 
ujz ^ i^±, i-e., cigar shaped. Following the successful approach of the NPSE description as for a single component 
scalar condensate we factor the wave function into transversal and longitudinal functions as 
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where x ^^'^ fi ^-re variational functions which depend on z and t. ^± is the transversal wave function, satisfying 
/ (irj_|$_Lp = 1, and is assumed identical for all components. Substituting Eq. |@J into Eq. (O, we obtain the 
Lagrangian of our system as 
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where V{z) = MujIz'^/2. E± is the transverse mode energy, and Ej_{x) = J df±^*^[-{h'^W']_/2M) + (i\/tj^r^/2)]$^. 
r] is the scaling factor of the nonlinear interaction strength, r]{x) = j dr±\'S?j_\'^, p{z) — and ^2 is independent 

of X and given by 
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To obtain the above result, we have also assumed a weak time and z dependence of the transverse wave function, i.e., 
d^±/dt ~ and V^$x — V^^^j^. The effective quasi-lD NPSE for a spin-1 condensate can now be derived from the 
least action principle of the above Lagrangian, 
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where p = pi is the total density and pi — is the density of the ith component. 



We discuss two separate ansatzes for the transverse 
function applicable respectively for the cases of weak and 
strong atomic interactions. 

• A Gaussian ansatz 

For weak atomic interaction when p — E± <^ l%uj± is 
satisfied, the transverse wave function can be taken 
as a Gaussian function of a variable width, 

$±(r±;x(^,t)) = ^^cxp[-ri/2x']. (7) 



The transverse mode energy and scaling factor are 
then given by 
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where a± = y^h/Muj±. 
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• A Thomas-Fermi ansatz 

For strong atomic interactions when fi~E± ^ huj± 
holds, the transverse wave function is taken as a 
Thomas- Fermi (TF) ansatz, 
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The kinetic energy in the transverse direction is 
neglected, leading to the transverse mode energy 
and scaling factor as 



4 



X 



(11) 
(12) 



CD 
Q 



0.1 

0.05 


0.04 
0.02 



N=10^ 


/' ^\ 


y 


/ / 
✓ y 




\\ 
\\ 
\v 

>^ 


-5 





5 


N=10^ 












-10 





10 



0.08 
0.06 
0.04 
0.02 



-% 



N=10' 




- - NPSE-TF 
■-■ NPSE-G 
3D 















10 



0.02 
0.01 




N=10^ 












-20 





20 



z/a 



z/a 



FIG. 1: The ground state density distribution of the conden- 
sate component in state |0) along the axis of the cigar-shaped 
trap, for *^Rb atoms and without an external magnetic field 
(flz = y^h/MuJz)- The inset shows the zoom-in central region. 
The solid line denotes the "exact" , while the dashed and dash- 
dot lines denote respectively the results from our NPSE with 
a TF or a Gaussian ansatz for the transverse profile. 

We performed some numerical simulations to illustrate 
the efficiency and effectiveness of the NPSE as developed 
by us for a spin-1 condensate in a cigar-shaped trap. For 
the first example, we computed the ground state of a 
^^Rb spin-1 condensate by propagating the GP equations 
and the effective ID NPSE with an imaginary time. The 
atomic parameters of ^''Rb are ag — 101.8 ob and 02 = 
100.4 Ob Ell. The trap frequencies are uj± = (27r)240 Hz 
and u!z = (27r)24 Hz. The "exact" solution as given by 
the ground state of the full 3D coupled GP equations © 
is calibrated by its effective ID distribution according to 
|/i(z)P — J df±\^i\'^. Figure n illustrates the results for 
several cases of different total number of atoms, TV. We 
note that with increasing N, the mean field interaction 
becomes stronger. For weak interactions the quasi-lD 
NPSE with a Gaussian variational ansatz gives a better 
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FIG. 2; The time dependence of the fractional condensate 
population in the |0) state Nq/N. The thick solid curve de- 
notes the full 3D simulation while the dashed curve denotes 
the simulation with our effective quasi-lD NPSE. We have 
used A'' — 10* for the top part and A'^ = 10^ for the bottom. 
As a comparison we also presented the result obtained from 
a time-independent Gaussian ansatz (thin solid curve in the 
bottom panel), which is shown to give a poor agreement in 
the strong interaction regime. 



result, while for strong interactions the quasi-lD NPSE 
with a TF ansatz is a better choice. Here "better" means 
the result obtained from an NPSE is closer to that of 
the full 3D solution. We also observe that the quasi- ID 
NPSE with a TF ansatz gives a lower central density 
and overestimates the TF radius in the weak interaction 
regime, while the quasi-lD NPSE with a Gaussian ansatz 
gives a lower central density and a correspondingly larger 
width in the strong interaction regime. Over all, it is in- 
teresting to point out that the quasi-lD NPSE with a 
Gaussian ansatz is not too bad even in the strong inter- 
action regime. 

To test the quasi-lD NPSE more strictly we study the 
dynamics of a spin-1 condensate out of equilibrium con- 
figuration and compare the results with those from a full 
3D simulation with the coupled GP equations ||2J). The 
initial state is taken as the ground state in a given mag- 
netic field. The simulation starts after the magnetic field 
is set zero, and we follow the spatial-temporal dynamics. 
With our NPSE, it becomes essentially a trivial task, 
and we find that excellent agreements are obtained with 
a Gaussian ansatz, e.g., with N — 10'^ atoms for weak in- 
teractions, and a TF ansatz for strong interactions with 
N = 10^ atoms. In the results to be given below, we 
instead use the effective quasi-lD NPSE to simulate the 
dynamics in the limit between the strong and weak inter- 
actions, i.e., for fj, — ?iu}± ~ h,uj±. For the quasi-lD NPSE 
approach, we use a Gaussian ansatz with N = 10^ atoms 
{£ - huj± = 1.74 S.w^) and a TF ansatz with TV = 10^ 
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FIG. 3: The contour plots of the density of the |0) component 
with respect to time and z. The left one is from the quasi-lD 
NPSE while the right one is the "exact" result from a full 3D 
simulation. Parameters are the same as Fig. |21 

atoms {S — Tlljo^ = 4.98 Tiwj,). Figure El displays the time 
evolution of the fractional condensate in the |0) state. For 
weak interactions, all three components share the space 
profile along the z axis, and the out of equilibrium dy- 
namics is periodic [iMllal23- FigureElalso clearly shows 
the periodic motion for N ~ 10^ atoms, although we do 
find that the quasi- ID NPSE gives a slightly shorter pe- 
riod than that of the full 3D simulation. For strong in- 
teractions, the apparent spatial profiles of the three spin 
components clearly become different, and the out of equi- 
librium dynamics also becomes complicated. Yet still, 
the quasi-lD NPSE simulations give results very close 
to the "exact" 3D solution, especially in the short time 
range. Figure |31 compares the dependence of the density 
distribution of the |0) state component on time and space 
from quasi- ID NPSE and full 3D simulation. The excel- 
lent agreement clearly demonstrates the efficiency and 
effectiveness of the quasi-lD NPSE approach, although 
we do find that it always seems to give a slightly shorter 
oscillation period as compared to the "exact" result. 



Before concluding, we hope to discuss the conditions 
under which our quasi-lD NPSE description is applica- 
ble. For a true ID condensate which enters the Tonks 
gas regime ,2L], our result is obviously not applicable. 
In the derivations of the quasi-lD NPSE, we have as- 
sumed a weak time- and z-dependence of the transverse 
mode. This assumption is valid only for weak excita- 
tions of the condensate such that the excitation in the 
transverse direction is negligible. In other words, the 
wavelength of the excitation is longer than the trans- 
verse size of the condensate. Under this condition, the 
transverse mode is reduced to the ground state which is 
a Gaussian in the weakly interacting limit and a TF pro- 
file in the strongly interacting limit. For spin-1 conden- 
sates widely discussed now, either of ^"^Na or ^^Rb atoms, 
the spin-dependent excitation is always weak since C2 
is two or three orders of magnitude smaller than cq. 
The wavelength of the spin wave is thus larger than the 
transversal size of the condensate for a cigar-shaped trap 
H H H H m 113. Our resulting NPSE model can thus 
be directly applied. For a strongly excited system, one 
has to include the transverse motion, with a more general 
approach as developed by Kamchatnov and Shchesnovich 
in |22| and Salasnich et al. in [23l |. 

In summary, we have extended the successful effective 
quasi-lD nonpolynomial Schrodinger equation (NPSE) 
for a single-component scalar condensate to a multi- 
component spin-1 condensate in a cigar-shaped trap. We 
have demonstrated its validity with a Gaussian ansatz 
for the transverse profile in the weak interaction regime 
and with a Thomas-Fermi (TF) ansatz in the strong 
interaction regime. We have further demonstrated its ef- 
fectiveness with studies on both the static (ground state) 
and dynamic properties of a spin-1 *^Rb condensate in a 
cigar-shape harmonic trap. With the effective quasi-lD 
NPSE, simulations for out of equilibrium condensate 
dynamics become rather efficient, thus allowing detailed 
comparisons with the recently observed spatial temporal 
dynamics. 
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